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Unstable spinor Bose-Einstein condensates are ideal candidates to create nonlinear three-mode interferome¬ 
ters. Our analysis goes beyond the standard SU(1,1) parametric approach and therefore provides the regime of 
parameters where suh-shot-noise sensitivities can be reached with respect to the input total average number of 
particles. Decoherence due to particle losses and finite detection efficiency are also considered. 
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Interferometers provide the most precise measurements in 
physics [1-3]. Hence, there is an urgent demand for novel the¬ 
oretical proposals and experimental techniques aimed at fur¬ 
ther increasing their sensitivity. Most of the current atomic 
and optical interferometers are made of linear devices such 
as beam splitters and phase shifters. Their phase uncer¬ 
tainty is fundamentally bounded by the shot-noise limit A0 ^ 
when using probe states made of average n uncor¬ 
related particles [4, 5]. It has been clarified that overcom¬ 
ing this bound requires engineering proper particle-entangled 
states [4] (see Refs. [6-8] for reviews). Using such states, 
sub-shot-noise (SSN) phase uncertainties have been demon¬ 
strated in several recent proof-of-principle experiments with 
atoms [9-14] and photons [15]. Yet, noise and decoherence 
limit the creation and use of quantum correlations [16]. It 
is therefore crucial to search for alternative schemes where 
probe states are classical and quantum correlations useful to 
reach SSN sensitivities are created inside the interferome¬ 
ter [9-11, 17]. 

In this Letter, we show that the coherent spin-mixing 
dynamics (SMD) in a spinor Bose-Einstein condensate 
(BEC) [18, 19] can be exploited to realize a nonlinear three¬ 
mode interferometer, as shown in Eig. 1. The SMD consists 
of binary collisions that coherently transfer correlated pairs of 
trapped atoms with opposite magnetic moment [22] from the 
rrif — 0 to the my = ±1 hyperfine modes, and vice versa. 
The probe state of the interferometer is classical, given by a 
condensate initially prepared in the my = 0 mode, and quan¬ 
tum correlations are created by the SMD. We hrst study the 
interferometer in the mean-held limit, the uif — 0 mode op¬ 
erator being replaced by a c-number. This analysis is valid for 
a large number of particles and low transfer rates. In this case, 
the interferometer operations belong to the SU(1,1) group and 
it is possible to obtain analytical predictions for the phase sen¬ 
sitivity. In optical systems, where transfer rates are rather 
low, the probe state needs to be very intense and the SU(1,1) 
approach is well justihed [20]. SU(1,1) optical interferome¬ 
try has been theoretically discussed [20, 23-25] and recently 
experimentally realized [26]. In contrast, experiments with 
spinor BECs [13, 27, 28] can be performed well outside the 
mean-held regime, with probe states of a relatively small num¬ 
ber of particles and - thanks to strong nonlinearities - compar¬ 


atively high transfer rates. We have thus also implemented a 
full three-mode quantum analysis. Within this framework, we 
can rigorously provide phase sensitivity bounds with respect 
to the average total number of particles h in input. Eor real¬ 
istic values of n, including particle losses and hnite detection 
efficiency, SSN is obtained in a regime where quantum cor¬ 
rections to the mean-held picture are important. 

Spin-mixing interferometry with BECs. The protocol out¬ 
lined in Eig. 1 follows hve steps: (I) probe state preparation 
- we consider empty my = ±1 modes and a BEC of average 
n atoms in the rrif — 0 mode, (II) a hrst SMD, (III) phase 
encoding, and (IV) a second SMD. Einally (V) the atoms are 
released from the trap: the three magnetic modes are spatially 
separated and the particle number is measured by imaging the 
atomic clouds. 

A standard description of the SMD is obtained in the single¬ 
mode approximation [29]: the condensate spatial wave func¬ 
tion ip{r) in the rrif = ±1 modes is assumed to be the same 
as in the rrif = 0 mode and it is given by the solution of 
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FIG. 1. Left: scheme of spin-mixing interferometry with spinor 
BEC, here represented in the / = 1 manifold. Right: when the 
rrif = 0 mode is treated parametrically (mean-held approach), 
the interferometer operations can be visualized on a hyperbolic 
surface by projecting the transformed state over SU(1,1) coherent 
states [20, 21]. 
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the Gross-Pitaevskii equation in the trapping potential [30]. 
This approximation is justified for a relatively low atom num¬ 
ber n < 10®, and tight confinement, when the spin healing 
length is larger than the size of the atomic cloud. These con¬ 
ditions are fulfilled in typical experimental setups [19]. The 
field operators are thus approximated by ^i(r) = ip{r)di, 
where di (aj) are annihilation (creation) operators for modes 
i = TO/ = 0,±1 obeying the boson commutation relations 
[di,d}j] = Sij (Ni = ajoi is the particle number operator). 
Up to terms proportional to the constant total particle number 
N = N_i -f TVq -f iV_|_i, the many-body Hamiltonian describ¬ 
ing the SMD in a dilute atomic cloud is [30] 

Hsmd = + e“^*'^aoaoa^ia^) -f 

+ xn{No-l){N+i + N_i). (1) 

The first term is identical to four-wave mixing in nonlinear op¬ 
tics [20, 31], where (j) is the relative phase between the mj = 
0 and TO/ = ±1 modes. The second term in Eq. (1) is a mean- 
field shift. The coupling x = — cq) / 

depends on the s-wave scattering lengths cq and C 2 of two 
bosons of mass M scattering in the total spin channels F = 0 
and F = 2, respectively [18, 32]. We indicate as 
[{xt,4>)2] the parameters for the first [second] SMD. Ex¬ 
perimentally, the SMD can be accurately controlled via mi¬ 
crowave dressing [19] and, in particular, switched off during 
phase acquisition. Neglecting interaction between particles 
during this stage, the (linear) phase shift Hamiltonian is 


Hps=hq{N+i+N_i), (2) 

where hq is the energy difference between the to/ = 0 and 
the TO/ = ±1 modes, see Fig. 1. The unitary transformation 
g-zilpstps/s encodes the phase shift 6 = 9+i -f 0_i = 2qtps^ 
where 9±i are the phases accumulated by the atoms in the 
TO/ = ±1 modes, relative to the ones in the to/ = 0 mode, 
during a time tps- For instance, the signal can be the second- 
order Zeeman shift due to a sufficiently strong magnetic field. 
Note that the first-order Zeeman shift, proportional to the net 
magnetization (equal to zero for our initial state), is conserved. 

The phase shift is estimated by measuring the number of 
particles in the to/ = ±1 modes at the end of the inter¬ 
ferometric sequence. We calculate the phase uncertainty as 
= 1 /\/toF( 0), the Cramer-Rao lower bound [6, 8, 
33], where to accounts for the repetition of independent mea¬ 
surements. 


F(o)= V 1 

d9 ) 


is the Fisher information (El) and P{N±i\9) is the condi¬ 
tional probability to measure N±i particles given the phase 
shift 9. A9cr is a saturable lower bound of phase uncer¬ 
tainty [6, 8, 33]. The FI can be experimentally extracted fol¬ 
lowing the method demonstrated in [12]. Alternatively, we 
can calculate the phase uncertainty from the error propagation. 


A6»ep = (^±i)out is the average num¬ 

ber of particles in output and (AiVd:i)out i^ '^he corresponding 
variance. This method is experimentally feasible but not al¬ 
ways optimal: we have A9ep/\/m > A^cr, in general [4, 8]. 

Mean-field approach. When the initial condensate con¬ 
tains a large number of particles and is weakly affected by 
the SMD, we can study the interferometer operations by re¬ 
placing do with y/fi. We introduce the operators = 
l{dl_id^_i + d+id-i), ky = Y^kpid^ii - d+ia_i), = 
|(d^p;^a+i -I- d^_id-i + 1), which belong to the SU(1,1) 
group and satisfy [kx,ky] = —ik^, [ky,kz\ = ik^ and 
[kz^kx:] = iky [20, 34]. Equations (1) and (2) thus be¬ 
come, up to a constant term, fTsMD = (2n — l)x^kz + 
2nxli{kxCOs2(j) -f kySm2(j>) and Hps = 2tiqkz, respec¬ 
tively. The interferometer protocol starts with vacuum in the 
TO/ = ±1 modes [Fig. 1(1)]. The first SMD g-'^^sMot/h 
[(xOi = Xty 4>i = 0] generates a Lorentz boost [35, 36] 
that amplifies the population in the to/ = ±1 modes 

where Af = (A+i -f iV_i)sMD [Fig. 1(11)]. The mean-field 
description is thus valid when [37] 

Xf —> 0, n —)• -foo, such that 0 < x^V^ "C 1. (5) 


The SMD generates a thermal distribution of perfectly cor¬ 
related atom pairs in the ±1 modes [22]: the two-mode 
squeezed-vacuum state [31], with variance (AA±i)g]yjj 3 = 
+ 1). The transformation rotates the state 

around the z axis of an angle 9 [Fig. l(III)]. The final oper¬ 
ation is a second SMD. This can be implemented either as 
an inverse Lorentz boost [i.e. (xi )2 = ~(x0i’ 

(j )2 = 0, as in Fig. 1(IV)], or by applying a 7r/2 phase shift to 
the TO/ = 0 mode followed by the transformation 
[i.e. (xf )2 = (x^)iz = 7’‘/2]. The latter is easier to be 
realized experimentally [38]. In both cases, the conditional 
probabilities are 


P{N±i\9) 


2[N'{Af + 2)(1 — 

[A/”(A/”-I- 2)(1 — COS0) -I- 2]''^±i+i ’ 


A direct calculation of Eq. (3) yields 


F{9) 


+ 2 ) 

N{N + 2) sin^ I -p 1 


cos^ 


(7) 


where N is given by Eq. (4). The FI reaches its maximum 
at 0 = 0. In this case, if (xi )2 = ~(x0i '^he two SMDs 
exactly compensate and the output to/ = ±1 modes are 
empty. Note also that (iV±i)out = Af{Af F 2) sin^ | and 
(A7V±i)2^t = iAN)l^jjSik f [(AiVj^j^Dsin^ 1 + 1]: error 
propagation saturates the Cramer-Rao lower bound, A0ep = 
A0CR- At 0 = 0, we obtain A0 cr = 1/^/mkfjATkk), 
which is below the shot noise, A0 cr < l/Vrnkf, calculated 
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FIG. 2. (a) The FI as a function of 9 (solid red line) and error propa¬ 
gation l/(A@)ep (dashed blue line). The shaded area is F{6) > n. 
Here n = 50 and rj = 0.2. (b) Fopt as a function of n. Dots 
are numerical results and solid lines are quadratic fits to the data: 
Fopt = a(? 7 )n^, for n 3> 1. Here r/ — 0.05 (red), 0.1 (blue), 0.15 
(green) and 0.2 (black). The shaded area is Fopt > n- The inset 
shows a{ri) as a function of 77 (dots). The solid line is a quadratic fit 
a{rj) ~ rj^ — l.Sl?;®; the dotted line is 0 ( 77 ) = if ■ 




considering only the average population in m/ = ±1 after the 
first SMD [20, 26]. We notice here that the shot noise should 
be calculated with respect to the total resources, i.e. the to¬ 
tal average number of particles n in the input state. However, 
such an analysis is impossible within the SU(1,1) framework. 

Full quantum approach. We have thus performed a full 
three-mode quantum analysis, investigating the regime of pa¬ 
rameters beyond Eq. (5). Thanks to the symmetry of the 
Hamiltonian (1), we can restrict ourselves to the Hilbert sub¬ 
space spanned by Fock states {| (Vg, = \k,M — 

2k,k)}, with 0 < yfc < [f J [39, 40]. We take p = 

Tfl— I®’ -^> (input) probe state. 

We numerically calculate F{9) for different values of the 
parameters n, rj and 9, where 77 = {N+i + fV_i)sMD/fi is 
the fraction of particles transferred from the to/ = 0 mode 
to the TO/ = ±1 modes after the first SMD. We mainly fo¬ 
cus on the case {xt )2 = ~(xi)i which, as shown below, is 
optimal. Overall, the FI as a function of 9 shows a behavior 
qualitatively similar to Eq. (7), with a maximum at 6 * = 0, see 
Fig. 2(a). A first important result is that, for proper values of 77 , 
the FI can be larger than n or, equivalently, A0 cr <1/ 'Jmn. 
In other words, it is possible to attain SSN uncertainties with 
respect to the average input number of particles. 

A scaling analysis of the FI as a function of n at the op¬ 
timal point 9 — 0 [we indicate Fopt = maxg F(9)] shows 
that Fopt ~ a(r])n^ asymptotically in n (in our simulations 
n < 1000), see Fig. 2(b). A fit gives a(r/) « 77^(1 — 1 . 877 ) 
in the case (xOs/ (xt)i = [see the inset of Fig. 2(b)]. We 
thus conclude that A9cr 1/n with a prefactor depending 
on 77 . 

Figure 3 is the main result of this Fetter. In panel (a) we 
show Fopt as a function of the ratio (x^) 2 /(x^)i’ for differ¬ 
ent values of 77 . For relatively large 77 , outside the mean- 
field regime, the curves are asymmetric around zero. The 
optimal interferometer configuration is reached for (xt )2 = 
— (xf)i, but SSN can be also obtained for positive values of 
(X'^)i/(xt)2' inverting the sign of x in the second SMD trans¬ 
formation, which might be experimentally difficult, is not nec- 


FIG. 3. (a) Fopt as a function of the ratio (xl) 2 /(x^)i- Solid lines 
refer to different values of 77 . Here n = 200. Panels (b) and (c) 
show the phase-sensitivity portrait in the ( 77 , n)-parameter space for 
(xf)i = —(x ^)2 and (xf)i = (xl) 2 , respectively. SSN phase un¬ 
certainties are obtained for n larger than a critical value ncriq) (dots, 
the solid line being a guide to the eye). The dashed red line in both 
panels is ncriv) = (1 — 2ri)/rf, obtained from a mean-field calcu¬ 
lation, which agrees with the numerics in the limit (5). In all panels, 
the shaded area indicates SSN. 


essary to reach SSN sensitivities. Figures 3(b) and 3(c) show 
the regime of parameters where SSN can be achieved, for 
(xf )2 = — (xf)i {xt )2 = (xf)i’ respectively. For fixed 

77 , a critical value n„{ri) exists such that A 0 cr < 
for n > ncr(??)- Deviations from the mean-field prediction, 
ncriv) = (1 ~ 2 ? 7 )/ 77 ^, can be appreciated for small fi, espe¬ 
cially for (xf) 2 /(xf)i > are relevant in current EEC 

experiments [13, 27, 28]. 

Particle loss and finite detection efficiency. According to 
Eq. (4), the SMD is unaffected by decoherence processes that 
happen on time scales much longer than ^ 1 /(xv/f)- For 
sufficiently large n and fast phase encoding, the nonlinear 
interferometer thus appears to be robust to one-body losses 
(relevant for the spin-mixing dynamics in the / = 1 mani¬ 
fold [41]). In fact, this dissipation source - due to inelastic 
collisions of the ultracold trapped atoms with the background 
thermal cloud, or by off-resonant light scattering in a dipole 
trap - has a density-independent rate. Conversely, recombina¬ 
tion losses - whose rate depends on fi - may strongly affect 
the interferometer sensitivity. We have thus simulated two- 
body losses in the to/ = 0 mode (relevant for the spin-mixing 
dynamics in the / = 2 manifold [13, 27]) using a Monte Carlo 
wave-function approach [42]. Fet 7 indicate the depletion rate 
during the SMD operation [i.e. {N^if)) = n/{l + 2ytn) for 
X = 0]. Figure 4(a) shows the regime of parameters ( 77 , fi) 
where SSN sensitivities can be found. The SSN region shrinks 
when increasing 7 /x and, in particular, no SSN is found for 
7 /x ^ 0.04. The branch structure of the SSN regions is ex¬ 
plained by the characteristic effects induced by particle losses 
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FIG. 4. (a) SSN region in {rj, n)-parameter space including two- 
body losses in the m/ = 0 mode, with loss parameter 7 /x = 0.01, 
0.02 and 0.03 (from outer to inner regions). The thick blue line, for 
7 = 0, is the same as the solid line in Fig. 3. (b) ?7 as a function of 
time, for different values of 7 /x (solid lines). Here n = 200 and the 
dashed line is Eq. (5). (c) Fopt as a function of n for 7 = 0.2 and 
different values of 7 /x. Dotted lines are guides to the eye. (d) SSN 
region in ( 7 , n)-parameter space obtained for detection noise a = 1, 
2, 5, 10 (from outer to inner regions). The thick blue line (cr = 0) 
corresponds to the solid line in Fig. 3. In panels (a) and (d) the dots 
are numerical results and the solid lines are guides to the eye. 


shown in Figs. 4(b) and 4(c). In Fig. 4(b) we plot ry as a func¬ 
tion of time, for different values of ■y/x- Losses decrease the 
transfer rate and place an upper bound to the achievable rj. In 
Fig. 4(c) we show the FI as a function of n. For n <C (x/27)^, 
the effect of losses can be neglected and we recover the scal¬ 
ing Fopt c>c {rjn}'^ of the noiseless case. For n > (x/ 27 )^ 
losses dominate and the sensitivity quickly degrades. For in¬ 
stance, in typical experiments with ®^Rb in the f — 2 mani¬ 
fold, the coupling strength is x ~ 0.5 Hz and we estimate a 
ratio 7/x ~ 10“^ — 10“^, well within our explored range. 

To model finite detection efficiency we consider a Gaus¬ 
sian convolution of the ideal output probabilities [43,44]. Re¬ 
sults for different values of the detection noise a are shown in 
Fig. 4(d). In typical experiments cr « 10, while a high detec¬ 
tion sensitivity cr « 1 has been discussed in Ref. [11]. In the 
regime (5) we can evaluate the FI from a convolution of proba¬ 
bilities ( 6 ). This allows for semianalytical calculations giving, 
to the leading order in l/ry and for cr > 1 , ncrljj) « 2<j/rf', 
which agrees with numerical calculations for n +00 and 
rj ^ 0. It predicts that ncr(p) shifts toward larger values 
when increasing cr, an expected behavior [25] that qualita¬ 
tively holds also outside the mean-field regime. 

Conclusions. We have studied a nonlinear three-mode inter¬ 
ferometer with spinor BECs. The nonlinear spin-mixing dy¬ 
namics not only splits the initial cloud but, differently from a 
linear beam splitter, it also creates, at the same time, quantum 
correlations among particles, necessary to overcome the shot- 
noise limit. Therefore, differently from linear interferometers, 


the nonlinear scheme discussed in this Letter can reach SSN 
phase uncertainties with classically correlated probe states. 
Accurate predictions of the phase sensitivity require a full 
three-mode quantum analysis, beyond the SU(1,1) (mean- 
field) approach. We have performed such an analysis and 
showed that it is possible to overcome the shot-noise limit 
with respect to the total average number of atoms in input. 
We also provide the regime of parameters where sub-shot- 
noise uncertainties can be achieved, including losses and fi¬ 
nite detection efficiencies. Our results pave the way to atomic 
ultrasensitive spin-mixing interferometry [38]. 
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